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Coherence properties of the down-converted beams generated in spontaneous parametric down- 
conversion are investigated in detail using an iCCD camera. Experimental results are compared with 
those from a theoretical model developed for pulsed pumping with a Gaussian transverse profile. 
The resuhs allow to tailor the shape of correlation area of the signal and idler photons using pump- 
field and crystal parameters. As an example, splitting of a correlation area caused by a two-peak 
pump-field spectrum is experimentally studied. 
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I. INTRODUCTION 

Light emitted from spontaneous parametric down- 
conversion in a nonlinear crystal is composed of photon 
pairs. Two photons comprising a photon pair are called 
a signal and an idler photon for historic reasons. The 
first theoretical investigation of this process has been 
done in year 1968 [H. Already this study has revealed 
that frequencies and emission directions of two photons 
in a pair are fully determined by the laws of energy and 
momentum conservations. For this reason, there occurs 
a strong correlation (entanglement) between properties 
of the signal and idler photons. In an ideal case of in- 
finitely long and wide nonlinear crystal and monochro- 
matic plane-wave pumping, a plane-wave signal photon 
at frequency cj^ belongs just to one plane-wave idler pho- 
ton at frequency uji that is determined by the conser- 
vation of energy. Emission angles of these photons are 
given by the momentum conservation that forms phase- 
matching conditions. Possible signal (and similarly idler) 
emission directions lie on a cone which axis coincides with 
the pump-beam direction of propagation. 

However, real experimental conditions have enforced 
the consideration of crystals with finite dimensions 0, Q , 
pump beams with nonzero divergence [3, as well as 
pulsed pumping . During this investigation, the ap- 
proximation based on a multidimensional Gaussian spec- 
tral two-photon amplitude has been found extraordinar- 
ily useful [U ■ The developed models have revealed 
that spatial characteristics of a pump beam are trans- 
ferred to certain extent to these of a photon pair gen- 
erated in a nonlinear crystal, especially in case of short 
crystals p"3 - [l5j . These models have also been useful in 
quantifying real effects in applied experimental setups 
utilizing photon pairs [igI . [TtI . They have also been re- 
cently extended to photonic [H, [l^ and wave-guiding 
[2Cl| - [25| structures. Also effects at nonlinear boundaries 
have been taken into account [1^, [13] ■ 

In this paper, we continue the previous investigations 
of spatial photon-pair properties [28r,'31j by experimental 
study of transverse profiles of the down-converted beams 
as well as correlation areas of the signal and idler pho- 



tons using an iCCD camera [32| - |34| . Special attention is 
paid to the role of pump-beam parameters. Experimen- 
tal results are compared with a theoretical model that 
considers Gaussian spectrum and elliptical pump-beam 
profile. We note that also sensitive CCD cameras have 
been found useful in investi gati ons of spatial properties 
of more intense twin beams |35l437l | . 

The paper is organized as follows. A theoretical model 
is presented in Sec. II. Sec. Ill brings theoretical anal- 
ysis of parameters of a correlation area as well as spec- 
tral properties of the down-converted fields. An experi- 
mental method based on the use of an iCCD camera is 
discussed in detail in Sec. IV. The experimentally ob- 
served dependence of parameters of the correlation area 
on pump-beam characteristics and crystal length is re- 
ported in Sec. V. Sec. VI is devoted to splitting of the 
correlation area and its experimental observation. Con- 
clusions are drawn in Sec. VII. 



II. THEORY 

The process of spontaneous parametric down- 
conversion is described by the following interaction 



Hamiltonian H\ 



int 



HnAt) = eoj drx^') :E(+)(r,t)Er^(r,<)Er^(r,i) 

V 

+ H.C., (1) 



where Ep^' is the positive-frequency part of the pump 
electric-field amplitude, whereas eI (e| ') stands for 
the negative-frequency part of the signal (idler) electric- 
field amplitude operator. Symbol x*-^-* means the second- 
order susceptibility tensor and : is shorthand for tensor 
reduction with respect to its three indices. Susceptibility 
of vacuum is denoted as Eq, interaction volume as V and 
H.c. substitutes the Hermitian-conjugated term. 

We further consider parametric down-conversion in a 
LilOs crystal with an optical axis perpendicular to the 
z axis of fields' propagation direction and type-I inter- 
action. The pump field is assumed to be polarized ver- 
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tically (it propagates as an extraordinary wave) whereas 
the signal and idler fields are polarized horizontally (they 
propagate as ordinary waves). In this specific configura- 
tion, scalar optical fields are sufficient for the description. 
The interacting optical fields can then be decomposed 
into monochromatic plane waves with frequencies LOa and 
wave vectors kai 

Ei+\r,t) = y"dk„£;(+)(k„)exp(ik,r-it^,<)+H.c.; 



a — p,s, I. 



(2) 



The signal and idler fields at a single-photon level 
have to be described quantally and so their spectral 
amplitudes e'^^ {ka) can be expressed as e'^^ {\ia) = 
isjhwal ^2£QcAna(y!a)aa{}f^o) using annihilation opera- 
tors aa(ka) that remove one photon from a plane- wave 
mode ko in field a. Symbol h stands for the reduced 
Planck constant, c is speed of light in vacuum, A trans- 
verse area of a beam, and Ua means index of refraction 
in field a. 

Under these conditions, the interaction Hamiltonian 
i?int in Eq. ([T]) takes the form p^ : 

^int(i) = y rfk.y dk, y dkp4+)(kp) 

X exp {i [w(kp) - w(ks) - w(ki)] t) 
X y" dr exp \-i (kp - k^ - k;) r] 



X a 



t(k,)al(k,)+H.c. 



(3) 



We have assumed in deriving Eq. ^ that the func- 
tion A„(ws,Wi) = -h^iX! s'jJiX^'^'^ I {IcA^ Usif^ s)ni(yii) is 
a slowly varying function of frequencies Ws and Wi and 
can be approximated by its value taken at the central 
frequencies and o;,^. 

A quantum state \^) of a generated photon pair can 
be obtained after solving the Schrodinger equation up to 
the first power of the interaction constant that results in 
the formula: 



dti/int (t)|vac) 



(4) 



|vac) means the vacuum state. Substitution of the in- 
teraction Hamiltonian iJjnt from Eq. ([3]) into Eq. ([4]) 
provides the following form for the quantum state 



I*) = y"rfk,y"dk,5(k„k,)ai(k,)a|(k,)|vac), 



(5) 



where the newly introduced two-photon amplitude S 
takes the form: 



5(k„k,) = a; J dkpEl+\kp)S{ujp - uJs - u;,) 

X / drexp [-i(kp - k^ - ki)r] (6) 
Jv 



and = — 27ri/?iA„(w°, aj°). We note that squared 
modulus |S'(ks,ki)p of the two-photon amplitude gives 
us the probability density of simultaneous generation of 
a signal photon with wave vector k^ and its twin with 
wave vector k^. 

Spectral resolution is usually not found in experiments 
with photon pairs and then the photon-pair coincidence- 
count rate is linearly proportional to the fourth-order cor- 
relation function Ggj defined as: 



sin(fs)sin(fi) 



X j duJ^UJ'l\h{LOs)h{LO,)\^\S{^,s,5s,^s,£.^,5^.0J^)?\ (7) 

S{^s-,5s,LOs,ii.,5i,u:i) = S'(ks,ki). The propagation di- 
rection of a photon is parameterized by radial emission 
angles (determining declination from the z axis) and 
azimuthal emission angles 5a (describing rotation around 
the z axis starting from the x axis); a = s,i (see also 
Fig. IS]). Functions hg and hi introduced in Eq. (O de- 
scribe amplitude spectral and/or geometrical filtering of 
photons in front of detectors. 

More detailed information is contained in intensity 
spectrum Ss of a signal field assuming photon pairs emit- 
ted into the fixed signal- and idler-photon directions 
given by angles <5s, 6, and 6^■. 



sin(^s)sin(^i)a;^|/i(ius)|' 



X / duJ^0Jf\h{uJ,)\^\S{^s,Ss,iOs,^^,S„UJ,)\'^. (8) 



If the signal-photon emission direction described by an- 
gles and Ss is not resolved, an integrated signal-field 
emission spectrum S'"* is observed: 



■it/2 



7r/2 



dSsSsiuJs;^s,Ss,^i,Si). 



(9) 

Similar formulas as given in Eqs. ^ and (|9]) can be de- 
rived also for the idler field. 

On the other hand excluding resolution in emission 
directions, spectral correlations between the signal and 
idler fields are characterized by a two-photon spectral 
amplitude $^.1 which squared modulus is defined as: 
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d5.. 



dis / d5, / di, 



sin(^s) s.\n{£,i)\h(ujs)h{u}i)\'^ 



(10) 



We further consider a Gaussian pump beam with the 
electric-field amplitude E^^ in the from: 

E^j^\v,t) = dujpAp{ujp)exp{ikp^z - iLOpt) 



Wpx{z_ 



■ exp 



exp 
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1 



exp 



Wpy{z) 

X exp[iCp(z)] 



exp 



-ih 



(11) 



as: 



|kp|. The functions Wpa, Rpa, and are defined 




Rpa{z) 



Cp{z) — [arctan(z/zp^) + arctan(z/Zpj^)]/2. (14) 

The function Ap introduced in Eq. (jlip gives the pump- 
field amplitude temporal spectrum. Constants Zp^ and 
zfj„ describe positions of waists with radii M^" and W^„ 
in the a; and y directions, respectively. Function Wpx(z) 
[Wpy{z)] gives a radius of the beam in the x [y] direction 
and with wavefront curvature Rpx{z) [Rpy{z)] at position 
z. 

We assume that the nonlinear crystal is sufficiently 
short so that changes of the pump- field amplitude E^p''^ in 
the transverse plane along the z axis can be neglected. In 
this case, the pump-field amplitude e'^^ can be charac- 
terized both by its temporal spectrum Apiujp) and spatial 
spectrum Fp{kpx, kpy) in the transverse plane: 



dkr, 



dkn 



£;^+)(r,<) = J d0JpAp{u!p) I u,r^px I u,r.py 

Fp(kpj., kpj^) exp(ikpj.a;) exp(zkpj^y) 
X exp(ikp2z) exp(— iWpi); 



(15) 



The spatial spectrum Fp corresponding to the Gaussian 
beam written in Eq. pip and propagating along the z 
axis can be expressed as: 

1 o 

Fp (kp2; , kpy ) 



Wpx{Z0)Wpy{Z0) WpxWpy 



X exp 



k2 

pa; 
* ' px 



exp 



k^ 

py 

py 



exp[<p(zo)], (16) 



where the position zq lies inside the crystal. Complex 
spectral half- widths Wpx and Wpy of the spatial spectrum 
in the transverse plane are given as follows; 



Wpa = 2. 



1 



T^p2,(zo) 2i?2Jzo)^ 



a = x,y. 



(17) 



In the following we consider a Gaussian chirped pump 
pulse which temporal amplitude spectrum Ap can be ex- 
pressed in the form: 



Apiujp) = ^p 



: exp 



4(1 zap) P 



(18) 



In Eq. (jl8l) . Tp denotes pump-pulse duration, Op stands 
for a chirp parameter, and ^p is the pump-field amplitude. 
We note that the amplitude width Awp (given as full 
width at 1/e of the maximum) of the pulse written in 
Eq. (dll) equals 4^(1 + al)/Tp. 

Considering the pump-field amplitude E^^ as given in 
Eq. P3|) the two-photon amplitude S defined in Eq. © 
can be recast into the form: 

5(^s,(5s,a;s,Ci,^j,'^j) = cAp{ujs+uJi) 
Fp{\^gx ~t" k^a:, ksy -(- kjy) 



X Lysine 



X exp < —I 



(19) 



sinc(a;) = sin(x)/x. In deriving Eq. (fT9|) . we have as- 
sumed that the crystal extents from z = —Lz to z = 0, 
Lz being the crystal length. The transverse profile of 
crystal is also assumed to be sufiiciently wide. 



III. CORRELATION AREA, SPECTRAL 
PROPERTIES 



Correlation area is defined by the profile of probabil- 
ity density of detecting a signal photon in the direction 
described by angles (^^ , 5s ) provided that its idler twin 
has been detected in a fixed direction given by angles 
(^i,(5i). In coherence theory, this probability is given 
by the fourth-order correlation function Gc, ^ defined in 
Eq. ([7]). Because the correlation function G^^i is usu- 
ally a smooth function of its arguments, it can be con- 
veniently parameterized using angular widths (given as 
full- widths at 1/e of maximum) in the radial (A^s) and 
azimuthal (AJ^) directions. In general, parameters of the 
correlation area depend on properties of crystal material 
as well as crystal length, pump-field spectral bandwidth, 
and transverse pump-beam profile. The last two param- 
eters allow to tailor characteristics of the correlation area 
in wide ranges. 

In the theoretical analysis of Sec. Ill, we use radial {£) 
and azimuthal (<5) angles inside a nonlinear crystal. The 
reason is that we want to exclude the effect of mixing 
in spatial and frequency domains at the output plane 
of the crystal in the discussion. However starting from 
Sec. IV radial (^) and azimuthal {5) angles outside the 
nonlinear crystal are naturally used in the presentation 
of experimental results. 

In radial direction, crystal length and pump-field spec- 
tral bandwidth as well as transverse pump-beam profile 
play a role. The dependence of radial width of the cor- 
relation area on the crystal length emerges through 
the phase matching condition in the z direction. This 
condition is mathematically described by the expression 
sinc(Ak^L^/2) in Eq. (HH); Ak^ = kp^-k^^-k^^. Actual 
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radial width is determined by this condition and conser- 
vation of energy (ujp = oJs+oJi). According to the formula 
in Eq. (fT9|) . the longer the crystal, the smaller the radial 
width. Analytical theory also predicts narrowing of the 
signal- and idler-field spectra with an increasing crystal 
length. If pulsed pumping is considered, the wider the 
pump-field spectrum, the greater the radial width and 
also the greater the signal-field spectral width (compare 
Figs. [TJ:, d with Figs [1^, b). This can be understood 
as follows: more pump-field frequencies are present in a 
wider pump-field spectrum and so more signal- and idler- 
field frequencies are allowed to obey the phase-matching 
conditions in the z direction and conservation of energy. 
In more detail and following the graphs in Figs. [It and d, 
signal-field photons with different wavelengths are emit- 
ted into different radial emission angles • Superposition 
of photon fields emitted into different radial emission an- 
gles then broadens the overall signal-field spectrum. 
It is important to note that all idler-field photons have 
nearly the same wavelengths which means that signal- 
field photons emitted into different radial emission angles 
use different wavelengths of the pulsed-pump spec- 
trum. The transverse pump-beam profile affects the ra- 
dial width through the phase-matching condition in the 
radial plane. This radial phase matching condition is an 
additional requirement that must be fulfilled by a gen- 
erated photon pair. Qualitatively, the more the pump 
beam is focused, the wider its spatial spectrum in radial 
direction and so the weaker the radial phase-matching 
condition. However, this dependence is quite small in 
radial angles, as follows from the comparison of graphs 
in Figs. [1^, b and Figs. [TJj, f. On the other hand, focus- 
ing of the pump beam leads to considerable broadening 
of the signal- and idler-field spectra in all radial emis- 
sion directions. Finally, if a focused pulsed pump beam 
is assumed (see Figs. [1^, h), broadening of the correla- 
tion area in radial direction as well as broadening of the 
overall signal- and idler-field spectra is observed due to 
a final pump-field spectral width. On the top, broad- 
ening of the signal- and idler-field spectra corresponding 
to any radial emission angle occurs as a consequence 
of pump-beam focusing. This behavior is related to the 
fact that indexes of refraction of the interacting fields are 
nearly constant inside the correlation area. We can say 
in general, that spectral widths of the signal and idler 
fields behave qualitatively in the same way as the radial 
width of correlation area. 

Comparison of the signal- and idler-field spectra in 
Figs. [TJ;, d valid for pulsed pumping with those in 
Figs. [T^, b for cw pumping leads to a remarkable ob- 
servation. Photon pairs generated into different signal- 
photon radial emission angles use different pump-field 
frequencies. There occurs spectral asymmetry between 
the signal and idler fields that originates in different de- 
tection angles considered; whereas the idler-field detec- 
tion angle is fixed, the angle of a signal-photon detection 
varies. This asymmetry determines the preferred direc- 
tion of the signal- and idler- field frequency correlations as 
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FIG. 1: Contour plots of signal- [S's(Ai,)] and idler-field 
[Si{Xi)] intensity spectra as they depend on radial signal- 
field emission angle ^s', idler-field emission angle is fixed. 
Spectra are determined for cw plane- wave pumping (a, b), 
pulsed plane- wave pumping (c, d, AAp — 2.8 nm), cw focused 
pumping (e, f, Wp'^ — 20 /im) and pulsed focused pumping 
(g, h, AAp = 2.8 nm, W°'f = = IVpV = 20 /xm) for 

Lz — 5 mm. 



they are visible in the shape of squared modulus |$s,iP 
of two-photon spectral amplitude introduced in Eq. pH)) 
[a large signal-field detector is assumed]. Contour plot 
of the squared modulus I'&is.iP of two-photon amplitude 
has a typical cigar shape. In cw case, the main axis of 
this cigar is rotated by 45 degrees counter-clockwise with 
respect to the Xi axis (see Fig. [2^) in order to describe 
perfect frequency anti-correlation. If pulsed pumping is 
taken into account, the cigar axis tends to rotates clock- 
wise; the broader the pump-field spectrum, the greater 
the rotation angle. Even states with positively correlated 
signal- and idler-field frequencies can be observed for suf- 
ficiently broad pump-field spectra (see Fig. ^) . We note 
that different dispersion properties at different propaga- 
tion angles have been fully exploited in the method of 
achromatic phase matching that allows to generate pho- 
ton pairs with an arbitrary orientation of the two-photon 
spectral amplitude [l^ \^ . 

The azimuthal width of correlation area is determined 
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FIG. 2: Contour plots of squared modulus I'l's.iP of two- 
photon spectral amplitude for a) cw and b) pulsed (AAp = 
8.5 nm) plane-wave pumping; = 5 mm. 



predominantly by a pump-beam transverse profile for ge- 
ometric reasons. To be more specific, it is the pump- 
beam spatial spectrum in azimuthal direction that affects 
the azimuthal extension of the correlation area through 
the phase- matching conditions in the azimuthal direc- 
tion. As material dispersion characteristics of the crystal 
are rotationally symmetric with respect to the z axis (sig- 
nal and idler fields propagate as ordinary waves) , the az- 
imuthal width of spatial pump-beam spectrum does not 
practically influence spectral properties of the signal and 
idler fields. 

We illustrate the dependence of correlation area on 
pump-beam focusing using a 5 mm long crystal and both 
cw and pulsed pumping in Fig. [3] We can see in Fig. [3^1 
that the signal-field azimuthal width ASs is inversely pro- 
portional to the width Wp'^ (full- width at 1/e of the 
maximum; Wp'^ = W^^ = Wp^^) of the pump-beam 
waist whereas the radial width A^s does not practically 
depend on the width Wp'^ of the pump-beam waist. This 
is caused by the fact that the phase-matching condition 
in the z direction is much stronger than that in radial 
direction for a 5-mm long crystal and so the radial width 
A^s is sensitive only to the pump-field spectral width 
in this case. Pulsed pumping gives a broader correla- 
tion area in radial direction as well as broader signal-field 
spectrum compared to cw case (see Fig. [Sja). Increasing 
pump-beam focusing releases phase-matching conditions 
and naturally leads to a broader signal-field spectrum. 

Contrary to the azimuthal width, the radial width A^^ 
depends on the pump-field spectral width AAp. The 
larger the pump-field spectral width AAp the greater 
the radial width A^^ and also the greater the signal- 
field spectral width Autg, as documented in Fig. 2] for 
a focused pump beam. We can also see in Fig. that 
the radial width A^s reaches a constant value for suffi- 
ciently narrow pump-field spectra. This value is deter- 
mined by the phase-matching condition in the z direction 
for the central pump-field frequency Up and so depends 
on the crystal length Lz (together with material disper- 
sion properties of the crystal) . The longer the crystal the 
smaller the radial width A^s . 

The above described dependencies allow to generate 
photon pairs with highly elliptic profiles of the correla- 
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FIG. 3: a) Radial (A^s, solid curves) and azimuthal (A5s, 
dashed curves) widths of correlation area and b) signal-field 
spectral width Aujs as they depend on width Wp'-^ of the 
pump-beam waist for pulsed (AAp — 5 nm, triangles) and cw 
(AAp = 0.03 nm, circles) pumping; Lz = 5 mm. Logarithmic 
scales on the x and y axes are used. 



tion area provided that the pump-beam profile in the 
transverse plane is highly elliptic. As an example, we con- 
sider a pump beam having Wp^-f" /Wp^-f" = 10. The depen- 
dence of the radial (A^g) and azimuthal (AJg) widths and 
signal-field spectral width Aw^ on the central azimuthal 
signal-photon emission angle Ssq is shown in Fig. [5] in 
this case. Whereas the radial and azimuthal widths are 
comparable for the azimuthal signal-field emission an- 
gle SsQ = 7r/2, their ratio ASs/A^s equals approx. 20 
for 5sQ — 0. Focusing the pump beam from 200 fim to 
20 /im in radial direction results in doubling the signal- 
field spectral width Aw^ as documented in Fig. [5]d (see 
also Figs. [T^ and e). 



IV. EXPERIMENTAL SETUP 

We have used a negative uniaxial crystal made of 
LilOa cut for non-critical phase matching, i.e. the op- 
tical axis was perpendicular to the pump-beam propaga- 
tion direction. We have considered crystals of two dif- 
ferent lengths {Lz—2 mm and 5 mm) pumped both by 
cw and pulsed lasers. As for cw pumping, a semicon- 
ductor laser Cube 405 (Coherent) delivered 31.6 mW at 
405 nm and with spectral bandwidth AAp — 1.7 nm. 
The second-harmonic field of an amplified femtosecond 
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FIG. 4: a) Radial width A^s of correlation area and b) signal- 
field spectral width Aujs as functions of pump-field spectral 
width AAp for a 5-mm (circles) and 10-cm (triangles) long 
crystal assuming a focused pump beam; W°-f = 200 ^m. 
Logarithmic scales on the x and y axes are used. 
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FIG. 5: a) Radial (A^s, solid curve) and azimuthal (A^s, 
dashed curve) widths of correlation area and b) signal-field 
spectral width Aws as they depend on central azimuthal 
signal-field emission angle Sso for a highly elliptic pulsed pump 
beam (W°/ = 200 iJ.m, W°/ = 20 urn, AAp = 5 nm); 
Lz — 5 mm. 



Ti:sapphire system (Mira+RegA, Coherent) providing 
pulses at 800 nm and ~250 fs long was used in the pulsed 
regime. The mean SHG power was 2.5 mW at the crystal 
input for a repetition rate of 11 kHz. Spectral bandwidth 
was adjusted between 4.8 and 7.4 nm by fine tuning of 
the SHG process. A dispersion prism was used to sep- 
arate the fundamental and SHG beams (for details, see 
Fig.©. 

Transverse profile of the pump beam and its divergence 
were controlled by changing the focus length of converg- 
ing lens LI or using a beam expander (BE2X, Thorlabs). 
The used focal lengths Jli of lens LI laid in the interval 
from 30 to 75 cm. As we wanted the pump beam to be 
as homogeneous as possible along the z axis, we chose 
the distance z^i between the lens LI and the nonlinear 
crystal such that the beam waist was placed far behind 
the crystal, i.e. z^i < Spatial spectrum of the pump 
beam in the transverse plane as a very important param- 
eter in our experiment was measured by a CCD camera 
(Lu085M, Lumenera) placed at the focal plane of a con- 
verging lens L3. Spatial spectra in horizontal and vertical 
directions were determined as marginal spectra and pa- 
rameters Wpx and Wpy characterizing their widths were 
found after fitting the experimental data. A fiber-optic 
spectrometer (HR4000CG-UV-NIR, Ocean Optics) was 
used to obtain the pump-beam temporal spectrum after 
propagation through the nonlinear crystal. 

The experiment was done with photon pairs degener- 
ate in frequencies {Xso — ^io — 800 nm) and emitted in 
opposite parts of a cone layer (the central radial emis- 
sion angle was 33.4 deg behind the crystal). As shown 
in Fig. [6}d the signal beam was captured directly by a 
detector whereas the idler beam propagated to the de- 
tector after being reflected on a high-reflectivity mir- 
ror. Both beams were detected on a photocathode of an 
iCCD camera with image intensificr (PI-MAX:512-HQ, 
Princeton Instruments). Before detection, both beams 
were transformed using a converging lens L2, one narrow- 
bandwidth and two high-pass edge filters. The geometry 
of the setup was chosen such that the lens L2 mapped 
the signal and idler photon emission angles to positions 
at the photocathode; the photocathode was placed in the 
focal plane of lens L2. For convenience, lenses L2 with 
different focal lengths {fL2— 12.5, 15, and 25 cm) were 
used. The applied bandwidth filter was 11 nm wide and 
centered at 800 nm. Edge fihers (Andover, ANDV7862) 
had high transmittances at 800 nm (98%) and blocked 
wavelengths below 666 nm. 

An active area of the photocathode in the form of a 
rectangular 12.36 mm wide (see Fig. [7]) was divided into 
512 X 512 pixels. Spatial resolution of the camera was 
38 /zm (FWHM) and its main limitation came from im- 
perfect contrast transfer in the image intensifier. In order 
to make data acquisition faster the resolution was further 
decreased by grouping 4 x 4 or 8 x 8 pixels into one super- 
pixel in the hardware of the camera. Consequently, sev- 
eral tens of camera frames were captured in one second. 
The overall quantum detection efficiency including com- 
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b) 

FIG. 6: Experimental setup used for the determination of 
angular widths: a) Entire setup that includes both cw and 
pulsed pumping as well as pump-beam diagnostics (for more 
details, see the text), b) Detail of the setup showing paths of 
the signal and idler beams. 



ponents between the crystal and photocathode was 7%, 
as derived from covariance of the signal and idler pho- 
ton numbers. Widths of the signal and idler strips are 
given by the bandwidth filter and lens L2 focal length. 
As for timing, a 10 ns long gate of the camera was used 
synchronously with laser pulses. In cw case, a 2 /is long 
gate was applied together with internal triggering. This 
timing together with appropriate pump-field intensities 
assured that the probability of detecting two photons in 
a single super-pixel was negligible. In other words, the 
number of detection events divided by quantum detec- 
tion efficiency had to be much lower than the number of 
super-pixels. 

Also the level of noise was monitored in the third nar- 
row strip; 1.82% of detection events came from noise. 
Detailed analysis has shown that 90 % of noise photons 
were red photons originating from fluorescence inside the 
crystal. Scattered pump photons contributed by 8.4 % 
and only 1.6% of noise counts were dark counts of the 



X s X i 




a) b) 



FIG. 7: Photocathode with registered photons after a) illumi- 
nation by light coming from 20 000 consecutive pump pulses, 
b) one pump pulse. The signal and idler strips image small 
sections of the cone layer and are slightly curved. The cur- 
vatures are oriented in the same sense in both strips because 
the idler beam is reflected on a mirror. 

iCCD camera. 

The experimental signal-idler correlation functions 
and gy in the transverse plane described by horizontal 
(x') and vertical {y') coordinates of the reference system 
in this plane have been determined after processing many 
experimental frames. The formula for the determination 
of correlation function g^ can be written as follows (see 
also Fig. Eb): 

p—l m— 1 l—l 

(20) 

In Eq. (I^U)) . p indexes frames {N gives the number of 
frames) and m (1) counts signal (idler) detection events 
[up to Mp (Lp) in the p-th frame]. Symbol x'^' {x'f) 
denotes horizontal position of the l-th detection in the 
signal (idler) strip of the p-th frame. Correlations in 
the vertical direction given by the correlation function 
gy can be determined similarly. The formula in Eq. (j20p 
takes into account all possible combinations of pairwise 
detection events. Only some of them correspond to de- 
tection of both photons from one pair. The remaining 
combinations are artificial in the sense that they do not 
correspond to detection of a photon pair. This poses the 
following restriction to the method. The number of ar- 
tificial combinations that occur at random positions has 
to be large enough in order to create a plateau in a 2D 
graph of correlation function gx{x'g,x[). Real detections 
of photon pairs are then visible on the top of this plateau 
(see Fig. [5] later). 

Cartesian coordinates and y'p j — s, i, in the trans- 
verse plane can be conveniently transformed into angles 
(3j and measured from the middle (a;^'^''"*, y^'^°"') of the 
jth strip and defined in Fig. [8] using the formulas: 

7, = arctan [{x'^ - "')//l2] , 

/3j = arctan [{y'^ - y^'^™')//L2 cos(7j)] , j = s, i; 
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FIG. 8: Sketch showing the geometry of signal and idler 
beams. The photon emission direction is described by ra- 
dial (^) and azimuthal (5) emission angles. In detector plane, 
cartesian coordinates x' and y' are useful. Photons propaga- 
tion directions are then conveniently parameterized by angles 
f5 and 7. 



(21) 

/l2 means the focal length of lens L2. Angles /3j and 
are related to radial and azimuthal angles and dj by 
the following transformation: 

Pj = arcsin [sin(^j) sin((5j)] , 

7j = arctan[tan(^j)cos((5j)] -,^-,-,dct, j = s,i, (22) 

where the radial angle ^j.dot describes the position of a 
detector in beam j. 



V. EXPERIMENTAL DETERMINATION OF 
PARAMETERS OF CORRELATION AREA 

In the experiment, spatial and temporal spectra of the 
pump beam have been characterized first. Typical re- 
sults are shown in Figs. Eli and b and have been used 
in the model for the determination of expected parame- 
ters of the correlation area. The correlation area, or more 
specifically its radial and angular profiles, have been char- 
acterized using histograms gx{x'g,x[) and gyiy'siUi)- His- 
togram gx{x'g,x'j) [gy{y's,y'i)\ gives the number of paired 
detections with a signal photon detected at position x'^ 
together with an idler photon registered at position x'.^ 
[y'j\ . These histograms usually contain experimental data 
from several hundreds of thousands of frames. As graphs 
in Figs, int and d show detections of correlated photon 
pairs lead to higher values in histograms gx and gy around 
diagonals going from upper-left to lower-right corners of 
the plots. Finite spreads of these diagonals have their ori- 
gin in non-perfect phase matching and can be character- 
ized by their widths Aa;^ and Ay^. Or more conveniently 
by uncertainties in the determination of angles (3s and 
7s; A/3s sa Ay'J fL2 and Ajs ~ Ax'JjL2- As detailed in- 
spection of the histogram {gy) in Fig.[3J; (d) has shown. 



cuts of this histogram along the lines with constant val- 
ues of x- (?/■) do not depend on the value of x- (?/■). This 
reflects the fact that idler photons detected at different 
positions inside the investigated area on the photocath- 
ode have identical (signal-photon) correlation areas. This 
allows us to combine the data obtained for idler pho- 
tons detected at different positions together and increase 
the measurement precision this way. This approach thus 
provides the radial cross-section {Gs,i)fi^ of the correla- 
tion area along the radial angle 7s as a mean value over 
all possible values of the signal-field azimuthal angle (3s. 
Moreover, consideration of different idler-photon detec- 
tion positions means averaging over the angles 7^ and 
(3i. The averaging is indicated by symbol (). Mathemat- 
ically, the radial cross-section {Gs,i)p^ expressed in the 
coordinate x'^ can be derived along the formula 

{Gs,)M) - E-9 K , (23) 

where the function x'^^'^{x'j} gives the central position 
(given as a locus) of the cut of the histogram g{x'g,x'^ for 
a fixed value of the coordinate x'^. In the theory, the ra- 
dial cross-section {Gs,i)i3^ is determined using the fourth- 
order correlation function Gs,i written in Eq. ([7]), substi- 
tution of angles ^s, i^s, Cji and by angles js,(3s,Ji, and 
Pi [inverse transformation to that in Eq. ([221)] and finally 
integration over the angles /3s, 7i, and (3i. Similarly, the 
azimuthal cross-section (Gs.i)^^ of the correlation area 
along the azimuthal angle (is arises after averaging over 
the angles 7s, ji, and (5i and can be determined by a 
formula analogous to that given in Eq. ([231) ■ The radial 
and azimuthal cross-sections {Gs,i)i3^ and {Gs,i)'y, corre- 
sponding to the pump beam with characteristics defined 
in Figs.in^ and b are plotted in Figs.jHt and f. Solid lines 
in Figs. ^ and f refer to the results of numerical model 
and are in a good agreement with the experimental data. 

The radial width (A7s);3^ (measured as full-width at 
1/e of the maximum) of radial cross-section {Gs,i)/3^ de- 
pends mainly on the pump-beam spectral width AAp. It 
holds that the greater the pump-beam width AAp the 
larger the radial width {Ajs)i3, as documented in Fig.fTUl 
for crystals 2- and 5-mm long. In the experiment, 11- 
nm wide frequency filters have been applied to cut noise. 
However, certain amount of photons comprising a pho- 
ton pair has also been blocked. According to the theo- 
retical model, this has also resulted in a small narrowing 
of the radial cut of the correlation area (compare solid 
and dashed curves in Fig. \W\i . The theoretical curve in 
Fig-dnihas been experimentally confirmed for several val- 
ues of the width Wp.]/ of pump-beam waist both for cw 
and pulsed pumping. 

On the other hand and in our geometry, it is the width 
Wp^f of the pump-beam waist that determines the angu- 
lar width {A(is)j^ of angular cross-section {Gs,i)j,- Pre- 
dictions of the model for 2- and 5-mm long crystals are 
shown in Fig. [TTl by a solid curve. This curve has been 
checked experimentally for several values of the width 
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FIG. 9: Typical measurement of a correlation area for pulsed 
pumping having 327,600 frames; = 5 mm. a) Spatial spec- 
trum of the pump beam determined in the focal plane of lens 
L3. b) Temporal intensity spectrum of the pump beam as de- 
termined by a spectrometer (diamonds), solid line represents 
a multi-peak Gaussian fit. c), d) Experimental histograms 
gx{x'g,x'i} (c) and gy{y's,v'i) (d). e), f) Experimental radial 
{{Gs,i)i3s, e) and azimuthal ((Gs,i)7s, f) cross-sections of the 
correlation area (rectangles, circles) together with theoretical 
predictions (solid lines). 



of the pump-beam waist both for cw and pulsed 
pumping. We note that these curves do not depend 
on the pump-beam spectral width AAp. We can see in 
Fig. [TT]that the measured points agree with the theoret- 
ical curve for smaller values of the width Wp^-^ . Larger 
values of the width W^^-^ lead to small angular widths 
(A/3s)^^ that could not be correctly measured because of 
the limited spatial resolution of the iCCD camera. 



15-, 
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FIG. 10: Radial width {A'ys)i3, as a function of pump-beam 
spectral width AAp. Experimental points have been obtained 
for Lz — 5 mm, Wpi^ > 140 fj,m (triangles), Lz — 5 mm, 
Wp/ < 140 fim (diamonds), Lz = 2 mm, W°/ > 140 fim 
(crosses), and Lz = 2 mm, Wpi/ < 140 (circles) both for 
cw and pulsed pumping. The theoretical model gives the same 
dependence for both crystal lengths Lz in cases without (solid 
curve) as well as with (dashed curve, 11-nm wide) spectral 
filters. 




O-h ■ , ■ , ■ 1 

400 800 1200 
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FIG. 11: Angular width (A/3s)7s as it depends on width W^y^ 
of the pump-beam waist for pulsed (crosses and triangles) 
and cw (diamonds) pumping. 2-mm (triangles and diamonds) 
and 5- mm (crosses) long crystals were used in experiment. 
The theoretical model predicts the same dependence for both 
crystals (solid curve). The dashed horizontal line indicates 
the resolution limit given by the camera super-pixel size. 



VI. ENGINEERING THE SHAPE OF A 
CORRELATION AREA 

As the above results have shown parameters of a corre- 
lation area can be efficiently controlled using pump-beam 
parameters, namely temporal spectrum and transverse 
profile. Even the shape of correlation area can be con- 
siderably modified. Splitting of the correlation area into 
two parts that occurs as a consequence of splitting of the 
pump-field temporal spectrum can serve as an example. 



10 



Using our femtosecond pump system, we were able to ex- 
perimentally confirm this behavior. We have generated a 
pump beam with the spatial spectrum given in Fig. [T2k . 
Its temporal spectrum containing two peaks as was ac- 
quired by a spectrometer is plotted in Fig. [T2b . The ex- 
perimental radial width (A7s)^^ of cross-section {Gsa)p^ 
given in Fig. [T2h shows that the two-peak structure of 
the pump-field spectrum resulted in splitting of the cor- 
relation area into two parts. On the other hand and 
in agreement with the theory, the angular cross-section 
{Gs.i)'ys iiot affected by the pump-field spectral split- 
ting (see Fig. [T^V For comparison, the theoretical pro- 
file of the correlation area given by the correlation func- 
tion Gs,i and appropriate for the pump-beam parameters 
given in Figs. [T^ and b is plotted in Fig. [T^ . It indi- 
cates a good agreement of the model with experimental 
data. Moreover, the squared modulus of theoret- 

ical two-photon spectral amplitude reveals that splitting 
of the correlation area is accompanied by splitting of the 
signal-field spectrum (see Fig. [TT[). 



camera. Single detection events in many experimental 
frames are processed and provide histograms from which 
cross-sections of the correlation area can be recovered. 
This method has been used for investigations of the de- 
pendence of parameters of the correlation area on pump- 
beam characteristics and crystal length. The experimen- 
tally obtained curves have been successfully compared 
with a theoretical model giving fourth-order correlation 
functions. Radial profile of the correlation area depends 
mainly on pump-field spectrum and crystal length. On 
the other hand, azimuthal profile of the correlation area 
is sensitive only to the transverse profile of the pump 
beam. Splitting of the correlation area caused by a two- 
peak structure of the pump-field spectrum has also been 
experimentally observed. 
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f) 



FIG. 12: Determination of a correlation area for pulsed pump- 
ing composed of two spectral peaks; Lz = 5 mm. a) Spatial 
spectrum of the pump beam, b) Temporal pump-field in- 
tensity spectrum (experimental points are indicated by dia- 
monds, solid line represents a multi-peak Gaussian fit), c) ,d) 
Experimental radial {{Gs,i)i3g, c) and azimuthal ({Gs,^)^^, d) 
cross-sections of the correlation area (rectangles and circles) 
together with theoretical predictions (solid line), e) Contour 
plot of the theoretical correlation function Gs,;. f) Contour 
plot of the squared modulus of the theoretical two- 

photon spectral amplitude. 



